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Abstract 

The long time behavior of a couple of interacting asymmetric ex- 
clusion processes of opposite velocities is investigated in one space di- 
mension. We do not allow two particles at the same site, and a colli- 
sion effect (exchange) takes place when particles of opposite velocities 
meet at neighboring sites. There are two conserved quantities, and 
the model admits hyperbolic (Euler) scaling; the hydrodynamic limit 
results in the classical Leroux system of conservation laws, even beyond 
the appearence of shocks. Actually, we prove convergence to the set of 
entropy solutions, the question of uniqueness is left open. To control 
rapid oscillations of Lax entropies via logarithmic Sobolev inequality 
estimates, the symmetric part of the process is speeded up in a suitable 
way, thus a slowly vanishing viscosity is obtained at the macroscopic 
level. Following the stochastic version of Tartar-Murat theory 

of compensated compactness is extended to two-component stochastic 
models. 

Key WORDS: hydrodynamic limit, hyperbolic scaling, systems of con- 
servation laws, compensated compactness 
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1 Introduction 

The main purpose of this paper is to derive a couple of Euler equations (hy- 
perbolic conservation laws) in a regime of shocks. While the case of smooth 
macroscopic solutions is quite well understood, see and [Tl], serious 
difficulties emerge when the existence of classical solutions breaks down. A 
general method to handle attractive systems has been elaborated in JH]) see 
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also 0] and [H] for further references. Hyperbolic models with two conser- 
vation laws, however, can not be attractive in the usual sense because the 
phase space is not ordered in a natural way. We have to extend some ad- 
vanced methods of PDE theory of hyperbolic conservation laws to stochastic 
(microscopic) systems. Lax entropy and compensated compactness are the 
main key words here, see ^Hj, [TT] . [T3] . [T3j . [271] . [2] for the first ideas, and 
the textbook [17] for a systematic treatment. The project has been initiated 
in [3], a full exposition of techniques in the case of a one-component asym- 
metric Ginzburg-Landau model is presented in [S]. Here we investigate the 
simplest possible, but nontrivial two-component lattice gas with collisions, 
further models are to be discussed in a forthcoming paper [H]. Since the 
underlying PDE theory is restricted to one space dimension, we also have 
to be satisfied with such models. The proof is based on a strict control of 
entropy pairs at the microscopic level as prescribed by P. Lax, L. Tartar 
and F. Murat for approximate solutions to hyperbolic conservation laws. A 
Lax entropy is macroscopically conserved along classical solutions, but the 
microscopic system can not have any extra conservation law, thus we are 
facing with rapidly oscillating quantities. These oscillations are to be con- 
trolled by means of logarithmic Sobolev inequality estimates, and effective 
bounds are obtainable only if the symmetric part of the microscopic evolu- 
tion is strong enough. That is why the microscopic viscosity of the model 
goes to infinity, i.e. the model is changed when we rescale it. Of course, the 
macroscopic viscosity vanishes in the limit and thus the effect of speeding 
up the symmetric part of the microscopic infinitesimal generator is not seen 
in the hydrodynamic limit. 

Unfortunately, compensated compactness yields only existence of weak 
solutions, the Lax entropy condition is not sufficient for weak uniqueness in 
the case of two component systems. That is why we can prove convergence 
of the conserved fields to the set of entropy solutions only, we do not know 
whether this set consists of a single trajectory specified by its initial data. 
Let us remark that [15] has the same difficulty concerning the derivation 
of the incompressible Navier-Stokes equation in 3 space dimensions. The 
Oleinik type conditions of weak uniqueness are out of reach of our methods 
because they require a one sided uniform Lipschitz continuity of the Riemann 
invariants of the macroscopic system, see £Q for most recent results of PDE 
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theory in this direction. It is certainly not easy to get such bounds at the 
microscopic level. 

The paper is organized as follows. The microscopic model and the macro- 
scopic equations are introduced in the next two sections. The main result 
and its conditions are formulated in Section 4. Proofs are presented in Sec- 
tion 5, while some technical details are postponed to the Appendix. 

2 Microscopic model 

2.1 State space, conserved quantities, infinitesimal generator 

We consider a pair of coupled asymmetric exclusion processes on the discrete 
torus, particles move with an average speed +1 and —1, respectively. Since 
we allow at most one particle per site, the individual state space consists 
of three elements. There is another effect in the interaction, something like 
a collision: if two particles of opposite velocities meet at neighboring sites, 
then they are also exchanged after some exponential holding times. We 
can associate velocities ±1 to particles according to their categories, thus 
particle number and momentum are the natural conserved quantities; the 
numbers of +1 and —1 particles could have been another choice. 

Throughout this paper we denote by T n the discrete torus Z/nZ, n £ N, 
and by T the continuous torus R/Z. The local spin space is S = {—1,0, 1}. 
The state space of the interacting particle system of size n is 

Q n := s Jn . 

Configurations will generally be denoted as 

uj_ := (a->j)jgT™ £ £l n , 

We need to separate the symmetric (reversible) part of the dynamics. This 
will be speeded up sufficiently in order to enhance convergence to local 
equilibrium also at a mesoscopic scale. The phenomenon of compensated 
compactness is materialized at this scale in the hydrodynamic limiting pro- 
cedure. So (somewhat artificially) we consider separately the asymmetric 
and symmetric parts of the rate functions r : S x S — > R+, respectively, 
s : S x S — > The dynamics of the system consists of elementary jumps 
exchanging nearest neighbor spins: (u>j, u>j+i) — ► (Uj,Uj +1 ) = (uj+i,ujj), 
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performed with rate Xr(uj, uJj+i) + Ks(ujj,Uj + i), where X,k > are fixed 
speed-up factors, depending on the size of the system in the limiting proce- 
dure. 

The rate functions are chosen as follows: 

r(l,-l) = 0, r(-l,l) = 2, 
r(0,-l)=0, r(-l,0) = l, 

r(l,0)=0, r(0,l) = l, 

that is the rate of collisions is twice as large as that of simple jumps, and 

r(uj,U3 J+1 ) = uj (1 - ujJ +1 ) + - wf) , 

where u>j~ := 1(^=1} , 0J~ := l{ w .=_i} and 1a denotes the indicator of a 
set A . The rates of the symmetric component are simply 

The rates r define a totally asymmetric dynamics, while the rates s define 
a symmetric one. The infinitesimal generators defined by these rates are: 

L n f{ui) ■.= k^ ) ^+i)(/(0^+^) - m) 

K n f(u>) = s(uj,«; j+1 )(f(e jt j +1 u) - /(a;)) , 

j6T n 



where Ojj is the spin-exchange operator, 



ojj if k = i 
Ui if k = j 
uj k if k^i,j. 



Recall that periodic boundary conditions are assumed in the definition of 
L n and K n . 

To get exactly the familiar Leroux system (J2J as the limit, the two con- 
served quantities, rj and £ should be chosen as 

rjj = i](ujj) := 1 — \u>j\ and £j = £(ujj) := Uj. 

The microscopic dynamics of the model has been defined so that Y2j £j an d 
r]j are conserved, we shall see that there is no room for other (indepen- 
dent) hidden conserved observables. In terms of the conservative quantities 
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we have 

r(uj,u j+1 ) = -(1 - r)j - + Vj+i + (1) 
The rate functions are so chosen that the product measures 

with one-dimensional marginals 

1 — p ± u 

7rp,u(0) = p, tt PiU (±1) = . 

are stationary in time. We shall call these Gibbs measures. The parameters 
take values from the set 

V:={(p,u) € [0,1] x [-1,1] : p+\u\ < 1}, 

and the uniform 7r n := 7r™/ 3 q will serve as a reference measure. Due to con- 
servations, the stationary measures 7r" are not ergodic. Expectation with 
respect to the measures 7r" will be denoted by E P)U (-). In particular, given 
a local observable Vi := v(ui- m , . . . ,Wi+ m ) with m fixed, its equilibrium 
expectation will be denoted as 

T(p,u) := E PtU (vi). 

The system of microscopic size n will be driven by the infinitesimal 
generator 

G n = nL n + n 2 aK n , 

where a = cr(n) is the macroscopic viscosity, the factor na(n) can be in- 
terpreted as the microscopic viscosity. A priori we require that <r(n) <C 1 
as n — > oo. A very important restriction, y/na(n) ^> 1 will be imposed on 
<r(n), see condition in subsection 14.21 

Let Pq be a probability distribution on fi n , which is the initial distribu- 
tion of the microscopic system of size n, and denote 

the distribution of the system at (macroscopic) time t. The Markov process 
on the state space O n driven by the infinitesimal generator G n , started with 
initial distribution pfi will be denoted by X™. 
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2.2 Fluxes 

Elementary computations show that the infinitesimal generators L n and K n 
act on the conserved quantities as follows, see (pQ). 

L n i]i= -if}(ui,u) i+1 ) =: -^i + ^i-i, 

L n ii= -<p(ui,u)i+i) +<j)(ui-i,Ui) =: + (pi-i, 

K n r]i= -i) s (uji,uj i+1 ) + ip s (u i - 1 ,uj i ) =:-'0| + ^|_ 1 , 

K n ii= -(j) s (u)i,u i+l ) + 4> s (oj i - 1 ,u)i) =: -(/>! + (f)f_ v 

where 

= r(w i ,c J ; i+ i) (77, - T)i+i) 

<j)j = r(ui,u i+ i)(£i-€i + i) 

= 2 + ~ 2 + ) + 2 ~ } + " &+ 1 ) ' 

Note that the microscopic fluxes of the conserved observables induced by 
the symmetric rates s(u>j,u>j + i) are (discrete) gradients of the corresponding 
conserved variables. 

It is easy to compute the macroscopic fluxes: 

V(p,u) := Ep^ipj) = pu 
$(p,u) := E P|1i (0j) = p + u 2 

3 Leroux's equation — a short survey 

Having the macroscopic fluxes © computed, the Euler equations of the 
system considered are expected to be 

dtp + d x (pu) =0 

d t u + d x ( P + u 2 ) =0. 

with given initial data 

u(0,x) = u (x), p(0,x) = p (x). (5) 
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This is exactly Leroux's equation well known in the PDE literature, see |17j . 
In the present section we shortly review the main facts about this PDE. The 
first striking fact is that such equations may have classical solutions only for 
some special initial data, in general shocks are developed in a finite time. 
Therefore solutions should be understood in a weak (distributional) sense, 
and there are many weak solutions for the same initial values. 

The following vectorial notations sometimes make our formulas more 
compact: 

i * r 




u:=\ , $ 



op ou 



( (p_ \ 

2 



dp 2 dpdu 

d 2 d 2 



\ dpdu du 2 ) 

We shall use alternatively, at convenience, the compact vectorial and the 
explicit notation. 

3.1 Lax entropy pairs 

In the case of classical solutions (jlj can be written as dtu + D(u)d x u = 0, 
where 



D{p,u) := V*(p,«) 



u p 
1 2u 



is the matrix of the linearized system. The eigenvalues of D are just 
A = \(p,u) := u + - {y/u 2 + 4p + u} , 
p = p(p,u) := u - - {^/u 2 + 4p - u} . 

This means that (J2J) is strictly hyperbolic in the domain 

{(p,u):p>0,u£R, (p,u) 7^(0,0)}, 

with marginal degeneracy (i.e. coincidence of the two characteristic speeds, 
A = p) at the point (p, u) = (0, 0). 

Lax entropy/flux pairs (S(u), F(u)) are solutions of the linear hyperbolic 
system VF(u) = V S(u)-V&(u) , that is dtS(u)+d x F(u) = along classical 
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solutions. This means that an entropy S is a conserved observable. In our 
particular case this reads 

[F^= P S' p + 2uS' u . 
or, written second order linear equation for S: 

pS'pp + uS" u - S" u = 0. (7) 

This equation is known to have many convex solutions, see JU]. We call an 
entropy/flux pair convex if the map (p, u) \— ► S(p, u) is convex. In particular, 
a globally convex Lax entropy/flux pair defined on the whole half plane 

R + x R is 
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2u 3 

S(p,u) := plogp+ y, := up + uplogp + — . 

Weak solutions of © are called generalized entropy/flux pairs. Rie- 
mann's method of solving second order linear hyperbolic PDEs in two vari- 
ables (see Chapter 4 of jSj) and compactness of T> imply that generalized 
entropy/flux pairs can be approximated pointwise by twice differentiable 
entropy /flux pairs. 

An entropy solution of the Cauchy problem Q, (JSJ) is a measurable 
function [0, T] x T B (t, x) i— > u(t, x) G M + x R which for any convex en- 
tropy/flux pair (S, F), and any nonnegative test function ip : [0,T]xT->I 
with support in [0, T)xT satisfies 
r-t f 

(dtip(t,x)S(u(t,x)) + d x <p(t,x)F(u(t,x))) dxdt 

+ / <p(0,x)S(u(p,x))dx > (8) 

JT 

Note that u) = ±p, F(p, u) = ipu, respectively, S(p, u) = ±u, F(p, u) = 
±(p + u 2 ) are entropy/flux pairs, thus entropy solutions are (a special class 
of) weak solutions. Entropy solutions of the Cauchy problem P|). form 
a (strongly) closed subset of the Lebesgue space L p ([0, T] x T, dt dx) =: L v tx 
for any p £ [1, oo). 

3.2 Young measures, measure valued entropy solutions 

A Young measure on ([0, T] x T) x T> is v = u(t, x; dv), where 

(1) for any (t, x) G [0, T] x T fixed, v(t, x\ dv) is a probability measure on 



JT 



8 



T>, and, 

(2) for any A C T> fixed the map (t, x) t— > x; j4) is measurable. 
Given a probability measure f on R + x R, we shall use the notation 



(u,f):= / f(v)u(dv). 
Jv 

The set of Young measures will be denoted by y. A sequence v n G y 
converges vaguely to v G 3^, denoted i/ n — i^, if for any / £ C([0, T] x T x D) 

Jim / (u n (t,x) , f(t,x,-)) dtdx = / (u(t,x) , f(t,x,-)) dtdx, 
'o Jt Jo Jt 



n— >oo 



or, equivalently, if for any test function 6 C([0,T] X T) and any g G C(P) 



n^oo 



lim / / (p(t,x){v n {t,x) , g) dtdx = I I ip(t,x)(v(t,x) , g) dtdx. 
'o jt jo 



The set y of Young measures will be endowed with the vague topology 
induced by this notion of convergence, y endowed with the vague topology 
is metrizable, separable and compact. We also consider (without explicitly 
denoting this) the Borel structure on y, induced by the vague topology. 

We say that the Young measure v(t, x; dv) is Dirac-type if there exists 
a measurable function u : [0, T] xT-»P such that for almost all (t, x) G 
[0, T] x T, v(t,x;dv) = 5 u ^ t ^(dv). We denote the subset of Dirac-type 
Young measures by U C y. It is a fact (see Chapter 9 of J7j) that 



y = CO (y) = co(W) = U, 

where 'co' stands for convex hull and closure is meant according to the vague 
topology. 

We say that the Young measure u(t, x;dv) is a measure valued entropy 
solution of the Cauchy problem (jl j) . iff for any convex entropy/flux pair 
(S, F) and any positive test function ip : [0, T] x T — > M + with support in 
[0,T) x T, 

/ / (dt(p(t,x){v(t,x) , S) + d x (p(t,x){v(t,x) , F)) dxdt 
Jo Jt 

+ / (p(0,x)(v(0,x), S)dx > (9) 
Jt 

holds true. Measure valued entropy solutions of the Cauchy problem Q), 
(jnj) form a (vaguely) closed subset of y. 
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Clearly, if u : [0, T] x T — > P is an entropy solution of the Cauchy 
problem |J3J), ©in the sense of (jHJ), then the Dirac-type Young measure 
v(t, x; dv) := S u (t,x)(dv) is a measure valued entropy solution in the sense of 
iPjl. The convergence of subsequences of approximate solutions to measure 
solutions is almost immediate by vague compactness, the crucial issue is to 
show the Dirac property of measure valued entropy solutions. This is the 
aim of the theory of compensated compactness. 

3.3 Tartar factorization 

A probability measure u(dp,du) on R 2 satisfies the Tartar factorization 
property with respect to a couple (Si,Fi), i = 1,2 of entropy/flux pairs 
if 

{u, S X F 2 - S 2 F 1 ) = {u, Si){v, F 2 ) - (u, S 2 ) (u, Fi> . (10) 

Dirac measures certainly posses this property, and in some cases, there is 
a converse statement, too. The following one-parameter families of en- 
tropy/flux pairs play an essential role in the forthcoming argument: 

S a (p, u) := p + au - a? , F a (p, u) := (a + u)S a (p, u) , 
S a (p, u) := |/9 + an - a?\ , F a (p, u) := (a + u)S a (p, u) , 

where the parameter, a £ M. . The case of (S a , F a ) is obvious because it is a 
linear function of the basic conserved observables and their fluxes. 

The pair (S a ,F a ) satisfies © in the generalized (weak) sense. This is 
due to the facts that the line of non-differentiability, p + au — a 2 = 0, is just 
a characteristic line of the PDE ©, and (S a ,F a ) coincides with (±S a ,±F a ) 
on the domains D± := {=t(p + au — a 2 ) > 0}. 

Lemma 1. Suppose that a compactly supported probability measure, v on 
M? satisfies (fT0|) for any two entropy /flux pairs of type (fTTj) , Then v is 
concentrated to a single point, i.e. it is a Dirac mass. 

Proof. This is Exercise 9.1 in J2J 5 where detailed instructions are also added. 
For Reader's convenience we reproduce the easy proof. Suppose first that 
S a = p + au — a 2 = v-a.s. for some a E R , then (v, p) + a(v, u) = a 2 ; 
let a\ and a 2 denote the roots of this equation. Since S ai (p, u) = implies 
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S a2 (p,u) = 0, u = a\ + 02 and p = —a\a2 z^-a.s. Therefore we may, and do 
assume that 

(v,F a ) (u,u\p + au-a 2 \) 

gia) := - — =— = a + — — ; ^7— 

(u,S a ) (u,\p + au-a 2 \) 

is well defined for all a G R . It is plain that g(a) — a is continuous, bounded, 

and g(a) — a — > (f, it) as a — > ±00. Applying (jl()j) to (S a ,F a ) and (S a ,F a ) 

we get 5(a) — a = (z/, u5 a ) / (z>, 5 a ) . On the other hand, from (fTU|) for (5 a , F a ) 

and (Sb, Fb) we get 

(b- a) ((u,S a S b ) - (v,S a ){v,S b )) = {u, S a )(u, uS b ) - (v, S b ){v, uS a ) , (12) 

thus dividing by (6 — a)(i/, S a ){v, Sb) , and letting b — > a we see that g is 
differ entiable, and g'(a) > 1, consequently g(a) = a + (z/, u) for all a £ M. 
This means that 

(u, pu) + a(u, u 2 ) — a 2 (v, u) = (v, p)(u, u) + a(u, u) 2 — a 2 (v, u) 

for all a £ M, whence z/(ii 2 ) = z/ 2 (u) . Substitute now u = (v,u) back into 
(fT2|) . Since b — a ^ (v,u) may be assumed, we have 

(z>, S a S 6 ) = (z/, S a )(u, S b ) , 

consequently (v,p 2 ) = (u, p) 2 . □ 

This lemma establishes that measure solutions satisfying Tartar's fac- 
torization property (|1()|) are, in fact, weak solutions. 

4 The hydro dynamic limit under Eulerian scaling 

4.1 Block averages 

We choose a mesoscopic block size I = l{n). A priori 

1 < l{n) < n, 

but more serious restrictions will be imposed, see condition |B|) in subsection 
14.21 and define the block averages of local observables in the following way: 
We fix once for ever a weight function a : M. — > R+. It is assumed that: 

(1) x i — > a(x) has support in the compact interval [—1, 1], 

(2) it has total weight J a(x) dx = 1, 
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(3) it is even: a(—x) = a(x), and 

(4) it is twice continuously differentiable. 

Given a local variable Vi its block average at macroscopic space x is 
denned as 

v n {x) = C»(a,x) := -^(^T 1 ) v r ( 13 ) 

3 

Note that, since I = l(n), we do not denote explicitly dependence of the 
block average on the mesoscopic block size I. 

We shall use the handy (but slightly abused) notation 

v n (t,x) :=v n (X t n ,x). 

This is the empirical block average process of the local observable v%. 

In accordance with the compact vectorial notation introduced at the 
beginning of Section El we shall denote 

and so on. 

Let £ (t, x) be the sequence of empirical block average processes of the 
conserved quantities, as defined above, regarded as elements of L\ x := 
L l ([0, T] X T). We denote by ¥ n the distribution of these in L\ x : 

F n (A) := P (£ n E A) , (14) 

where A € L\ x is (strongly) measurable. Tightness and weak convergence 
of the sequence of probability measures ¥ n will be meant according to the 
norm (strong) topology of L\ x . Weak convergence of a subsequence P n ' will 
be denoted F n ' P. 

Further on, we denote by v n the sequence of Dirac-type random Young 
measures concentrated on the trajectories of the empirical averages £ (t, x) 
and by Q n their distributions on 3^ 

u n (t,x;dv) := S r{t ^dv), Q n (A) := P {v n E A) , (15) 

where A £ y is (vaguely) measurable. Due to vague compactness of y, the 
sequence of probability measures Q n is automatically tight. Weak conver- 
gence of a subsequence Q™' will be meant according to the vague topology 
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of y and will be denoted Q n =^ Q. In this case we shall also say that 
the subsequence of random Young measures v n> (distributed according to 
Q n ) converges vaguely in distribution to the random Young measure v (dis- 
tributed according to Q), also denoted v n =^ v. 

4.2 Main result 

All results are valid under the following conditions 

(A) The macroscopic viscosity a = o~(n) satisfies 

n~ 1/2 < a < 1. 

(B) The mesoscopic block size I = l(n) is chosen so that 

n 2/3 a 1/3 « I < na 



(C) The initial density profiles converge weakly in probability (or, equiv- 
alently in any L p , 1 < p < oo). That is: for any test function 
tp : T -> M x R 



lim E 

n— >oo 



<£>(x) • (^™(0,x) - u Q (x)) dx 



0. 



Our main result is the following 

Theorem 1. Conditions PI) . fBl) . and f£j) are m force. The sequence of 
probability measures P n on L\ x , defined in is tight (according to the 
norm topology of L\ x ). Moreover, if¥ n ' is a subsequence which converges 
weakly (according to the norm topology of L\ x ), P n => P, then the limit 
probability measure P is concentrated on the entropy solutions of the Cauchy 
problem ^jj, |3J). 

Remark: Assuming uniqueness of the entropy solution u(t, x) of the Cauchy 
problem (J3J), Q, we could conclude that 

£ — u, in probability. 



13 



5 Proof 



5.1 Outline of proof 

We broke up the proof into several subsections according to what we think 
to be a logical and transparent structure. 

In subsection 15.21 we state the precise quantitative form of the conver- 
gence to local equilibrium: the logarithmic Sobolev ineqaulity valid for our 
model and Varadhan's large deviation bound on space-time averages of block 
variables. As main consequence of these we obtain our a priori estimates: 
the so-called one-block estimate and a version of the so-called two-block es- 
timate, formulated for spatial derivatives of the empirical block averages. 
These estimates are of course the main probabilistic ingredients of the fur- 
ther arguments. The proof of these estimates is postponed to the Appendix 
of the paper. 

In subsection 15.31 we write down an identity which turns out to be the 
stochastic approximation of the PDE Various error terms are defined 
here which will be estimated in the forthcoming subsections. 

In subsection 15.41 we introduce the relevant Sobolev norms and by using 
the previously proved a priori estimates we prove the necessary upper bounds 
on the apropriate Sobolev norms of the error terms. 

In subsection 15.51 we show that choosing a subsequence of the random 
Young measures ()15JI which converges vaguely in distribution, the limit (ran- 
dom) Young measure is almost surely measure valued entropy solution of 
the Cauchy problem Q, ©• 

Subsection 15.61 contains the stochastic version of the method of com- 
pensated compactness. It is further broken up into two sub-subsections as 
follows. In sub-subsection IB . 6 . 1 1 we preent the stochastic version of Murat's 
Lemma: we prove that for any smooth Lax entropy/flux pair the entropy 
production process is tight in the Sobolev space Hf£- in sub-subsection 
15.6.21 we apply (an almost sure version of) Tartar's Div-Curl Lemma lead- 
ing to the desired almost sure factorization property of the limiting random 
Young measures. Finally, as main consequence of Tartar's Lemma, we con- 
clude that choosing any subsequence of the random Young measures (fT3)l 
which converges vaguely in distribution, the limit (random) Young measure 
is almost surely of Dirac type. 

The results of subsection 15.51 and sub-subsection 15.6.21 imply the Theo- 
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rem. The concluding steps are presented in subsection 15.71 



5.2 Local equilibrium and a priori bounds 

The hydrodynamic limit relies on macroscopically fast convergence to (local) 
equilibrium in blocks of mesoscopic size I. Fix the block size I and (N, Z) £ 
N x Z with the restriction N + \ Z\ < I and denote 

I l 

n l NiZ -.= {uien 1 : j> = w, j> = z}, 

3=1 3=1 

I I 

3=1 3=1 

and, for / : £l l N z — > E 

== E(/( w+^)-/^))> 



3=1 

-i 



3=1 

In plain words: 17^ ^ is the hyperplane of configurations w 6 SI 1 with fixed 
values of the conserved quantities, tt 1 n z is the microcanonical distribution 
on this hyperplane, K l N z is the symmetric infinitesimal generator restricted 
to the hyperplane ^^v^, and finally D l N z is the Dirichlet form associated 
to K l N z . Note, that K l N z is defined with free boundary conditions. Expec- 
tations with respect to the measures 71^ z are denoted by z ( • ) . The 
convergence to local equilibrium is quantitatively controlled by the following 
uniform logarithmic Sobolev estimate: 

Lemma 2. There exists a finite constant H such that for any I £ N, (N, Z) G 
N x Z with N + \Z\ < I and any h : n l NZ -> M + with E l N z (h) = 1 tfie 
following bound holds: 

E l NjZ (h log fc) < « Z 2 !^ (Vh) . (16) 



Remark: In (see also J2|) the similar statement is proved (inter alia) 
for symmetric simple exclusion process. That proof can be easily adapted to 
our case. Instead of stirring configurations of two colours we have stirring 
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of configurations of three colours. No really new ideas are involved. For 
sake of completeness however, we sketch the proof in subsection 16.11 of the 
Appendix. 

The following large deviation bound goes back to Varadhan . See 
also the monographs (Oj and @]- 

Lemma 3. Let I < n, V : S — > M+ and denote Vj(u£) := V(u>j, . . . ,Uj+i-i). 
Then for any (3 > 



Remarks: (1) Assuming only uniform bound of order l~ 2 on the spectral 
gap of K l N z (rather than the stronger logarithmic Sobolev inequality Ijlfijl) 
and using Rayleigh-Schrodinger perturbation (see Appendix 3 of we 
would get 



which wouldn't be sufficient for our needs. 

(2) The proof of the bound (|TT|) explicitly relies on the logarithmic Sobolev 
inequality p6|). It appears in j^H] and it is reproduced in several places, see 
e.g. [UE1- We do not repeat it here. 

The main probabilistic ingredients of our proof are the following two 
consequences of Lemma |3J These are variants of the celebrated one block 
estimate, respectively, two blocks estimate of Varadhan and co-authors. 

Proposition 1. Assume conditions LT]) and fS)) . Given a local variable 
Vj there exists a constant C (depending only on Vj) such that the following 
bounds hold: 




(17) 






(18) 




Ca 



(19) 
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The proof of Proposition Q is postponed to subsection 16,31 in the Ap- 
pendix. It relies on the large deviation bound (|17|) and an elementary prob- 
ability lemma stated in subsection 16.21 of the Appendix. 

We shall refer to (|18|) as the block replacement bound and to (|19|) as the 
gradient bound. 

5.3 The basic identity 

Given a smooth function / : D — > K we write 

dtf(T(t,x)) = G n f(T(t,x)) + d t Mf(t,x), 

where the process t \— > Mj(t,x) is a martingale. Here and in the future 
dtf(£ (t, x)) and dtM^(t, x) are meant as distributions in their time variable. 

In this order we compute the action of the infinitesimal generator G n = 
nL n + n 2 aK n on /(£ (x)). First we compute the asymmetric part: 

nL"/(f(*)) = -V/(r(^))-5,0 n (rE)+4 n (rE) (20) 

where 

A 1 j ,n (x) = A^' n (uj_,x) := n^^r(u>j,ujj + i) x (21) 

1ST 

{/(Too - yK^r 1 ) - °( " g 7" 1 ) ) " ~ 

+ ^'(^)v/(? n w)-fe-^, +1 )}. 

See formula © for the definition of cb. A l J n is a numerical error term which 
will be easy to estimate. 

Next, the symmetric part: 

n 2 aK n f(T(x)) = aVf(T(x)).d 2 x 'i n (x)+A 2 / n (x) (22) 

where 

A 2 f ' n (x) = A 2 f(u,x) := (23) 
-/(?>)) +^"(^)V/(?>))-^}. 
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This is another numerical error term easy to estimate. 
Hence our basic identity 



-^-7% -^-TL '^■Ti -^Tt 

dtttt (t, x)) + V/(£ (t, x)) • V<&(£ (t, x)) ■ d x £ (t, x) = (24) 



B) n {x) = B} n {ui,x) := d x {v/(f(x)) • {*(T(x)) - $>))} (25) 

Bf\x) = B} n (u,x) := ad 2 J(T(x)) = d x {<xV/(f (*)) ■ d x T(x)} (26) 

Cf{x) = C)> n {u>,x) := -{d x T(x)) ] ■ V 2 /(f(*)) • (*(?>)) - 4> n {x)) 

(27) 

Cf n (x) = Cf n (u,x) := -a{d x T(x)y ■ V 2 /(f (*)) ■ (d x T(x)) (28) 



In the present paper we shall apply the basic identity ()24|) only for Lax 
entropies f(u) = S(u). In this special case the left hand side gets the form 
of a conservation law: 



2 



i=l 



The various terms on the right hand side are 



and 



Af{t,x) :=Af(Xp,x), 
B l f(t,x) :=Bf{X^x\ 
Cfj n (t,x) :=C) n (X^x). 



d t S(T(t,x)) + d x F(T(t,x)) = 



(29) 



2 



]T (Af(t,x) + B^(t,x)+C^(t,x))+d t M'g(t,x), 



5.4 Bounds 



We fix T < oo and use the LP norms 




and the Sobolev norms 



:=sup{ / ip(t,x)g(t,x)dxdt : \\d t ip\\ 9 Tq + \\d x (p\\ q q < l} 
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where p 1 + q 1 = 1 and ip : [0, T] x T — > M is a test function. We use the 
standard notation W tx ' =: H tx . 

Remark on notation: The numerical error terms A l j n (t, x), i = 1,2, will 
be estimated in L^° x norm. In these estimates only Taylor expansion bounds 
are used, no probabilistic argument is involved. The more sophisticated 
terms B^ n (t,x), i = 1,2, respectively, C^ n (t,x), i = 1,2, will be estimated 
in H^ x , respectively, x norms. The martingale derivative dtMJ(t, x) will 
be estimated in H^ x norm. 

By straightforward numerical estimates (which do not rely on any prob- 
abilistic arguments) we obtain 

Lemma 4. Assume conditions iLTj) and W\) . Let f : T> — > M. be a twice 
continuously differentiable function with bounded derivatives. Then almost 
surely 



A 



1,7! 
/ 



o(l) and 



,2,n 



0(1) 



as n — > oo. 



Proof. Indeed, using nothing more than Taylor expansion and boundedness 
of the local variables we readily obtain 



sup sup 

xgt wen™ 



sup sup 

xeTuien n 



A) n {ui,x) <C~ = o(l) 



A 2 f ' n (ui,x) 



< C 



n 2 a 



o(l). 



(30) 



(31) 



We omit the tedious but otherwise straightforward details. 



□ 



Applying Proposition^we obtain the following more sophisticated bounds 



Lemma 5. Assume conditions iLTj) and W\) . Let f : T> — > M. be a twice 
continuously differentiable function with bounded derivatives. The following 
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asymptotics hold, as n — > oo: 



l,n 



2,» 



l.n 



2,d 



h; 



H t,x 



L t,x 



L t,x 



(i) E( B 

(«) 
(tit) 

(iv) 
Proof. 

(i) We use the block replacement bound Q18J) : 

r-T 

'0 Jt 



-o(l) 
-o(l) 
o(l) 
0(1) 



E 




v{t,x)B 1 j' n (t,x) dx dt ^ 



E 








d x v(t,x)Vf(T(t,x)) • (*(£"(*, *)) - 0* (*,*)) dsd* 



r-T 




< sup |V/(u)| ||^|| L? E 
ixgx> ,j: k Jo Jt 

<c||M| i? 

(ii) We use the gradient bound ()19j) : 



*(£ (t,x))-<f> (t,x) 



dx dt 



1/2 



E 




v(t, x)B'j ,n (t, x) dx dt 



to Jt 
E 



cT 




o Jt 



d x v(t,x)Vf(£ n (t,x)) ■ a(d x i n (t,x)) dxdt 



< sup |V/(«)| \\d x v\\ L f o-B 

uev ' x x 7o Jt 



r-T 




da; (it 



1/2 



<C\\d x v\\ L , a 1 ' 2 . 

t,x 

(Hi) We use both, the block replacement bound (j!8j) and the gradient bound 
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m 

E 



T 

JT 



dx dt 



<sup|V 2 /(«)|E( f T [ 4> n (s,x)-*(r(s,x)) 



2 x 1/2 

dxdt) x 



E 



T 



o Jt 



d x £ (s,x) 



dx dt 



1/2 



< c- 



na 



(iv) We use again the gradient bound (fl9|) : 



E 



o Jt 



Cf n (t,x) 



dx dt 



< sup |V 2 /(w)| o-E 
uev 

< C. 



JT 



dxi (s,x) 



dx dt 



□ 

Lemma 6. Assume conditions and Let f : T> — > M. be a twice 

continuously differentiate function with bounded derivatives. There exists a 
constant C ( depending only on f) such that the folowing asymptotics holds 
as n — > oo: 



Proof. Since 



E (ll^ M /L-)=°( 1 ) 



\d t Mf\\ 2 , < \\MJ\\% , 



we have to bound the expectation of the right hand side. 



E 



(J j (Mf(t,x)) 2 dxdt} = B (J o j ' {Mf(t,x))dxdtj, 



where t \— > (MJ(t,x)) is the conditional variance process of the martingale 
Mf(t,x): 



(M?(t, x)) = n [L n f\d (t, x)) - 2/(| (t, x))L n f{£ (t, x)) 

+n 2 a (K n f\T(t,x)) - 2f(T(t,x))K n f(T(t,x))) . 
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Using the expressions 1)20(1 and 1)22(1 we obtain 

(Mf(t,x)) = A 1 f(t,x)-2f(T(t,x))A 1 f > n (t,x) 
+A 2 f(t,x)-2f(i n (t,x))A 2 f(t,x). 

Hence, by the bounds ((HO)) and l)3T)l (which apply as well of course to the 
function f 2 ), we obtain 



.n 2 a 



sup sup(M?(t,x)) < C-==- = o(l), 

tG[0,T] ic€T ' 

which proves the lemma. □ 

5.5 Convergence to measure valued entropy solutions 

Proposition 2. Conditions Oj) . fBj) . and f7^) are in force. Let Q n be a 
subsequence of the probability distributions defined in U5\) . which converges 
weakly in the vague sense: Q™' =^ Q. Then the probability measure Q is 
concentrated on the measure valued entropy solutions of the Cauchy problem 

Proof. Due to separability of C([0,T] x T) it is sufficient to prove that for 
any convex Lax entropy/flux pair (S,F) and any positive test function tp : 
[0,T] xT-t R + , © holds Q-almost-surely. So we fix (S,F) and 93, and 
denote the real random variable 

X n := -/ hp(t,x)(d t S(T(t,x)) + d x F^ (t,x)))dxdt 
Jo Jt 

T r 

(dt<p(t, x) (v n (t, x) , S) + d x tp(t, x) (v n (t , x) , F)) dx dt 



JT 



+ 



[ <p(0,x)(v n (0,x) , S) dx. 
Jt 



The right hand side is a continuous function of v n , so from the asumption 
Q n =\ Q it follows that 

X n X, (32) 

where 

X := jj (d t (p(t,x)(u(t,x) , S) + d x ip(t,x)(v(t,x) , F)) dxdt 
Jo Jt 

+ / / v(0,x;dv)p(0,x)S(v)dx. 
Jt Jv 
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and v is distributed according to Q. 

We apply the basic identity pi]) speciafied for f(u) = S(u), that is 
identity ((21 • It follows that 

X n = Y n + Z n (33) 

where 

Y n := I [ ^(t,x)C 2 s ' n (t,x)dxdt 
Jo Jt 

= a ( T [ <p(t,x)(d x T(t,x)y -V 2 S(T(t,x))- (dJT&x)) 
Jo Jt 

and 

Z n ■= £ J <p(t, x) ( A T + B s n ) + c s n + Ws^j (t, x) dx dt. 

Due to convexity of S and positivity of (p we have 

Y n > 0, almost surely. (34) 

On the other hand, from Lemmas 0J El El we conclude that 

lim E(|Z n |) =0. (35) 

Finally, from (|3*2*|) . (|3*3~)) . and (|33|) the statement of the Proposition 
follows. □ 

5.6 Compensated compactness 
5.6.1 Murat's lemma 

Lemma 7. Assume conditions QJ) and Given a twice continuously 

differentiable Lax entropy / flux pair (S,F), the sequence 

X n (t, x) := d t S(i n (t, x)) + d x F(i n (t, x)) 
is tight in . 

Proof. Note that X n (t,x) is exactly the left hand side of the basic identity 
(j29j) and racall that this expression (in particular dtS(£ (t, x))) is a random 
distribution in its t variable. 
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By definition and a priori boundedness of the domain T>, there exists a 
constant C < oo such that 

P(\\x n \\ w -^<C) = 1. (36) 

We decompose 

X n (t,x) = Y n (t,x) + Z n (t,x), (37) 

where 

Y n (t, x) := B l s > n (t, x) + B% n {t, x) + d t M%(t, x), 

Z n {t, x) := A^ n (t, x) + A 2 ^ n (t, x) + Cl' n (t, x) + C 2 s ' n (t, x). 

For the definitions of the terms A l £ n , B^ n , Cg™, % = 1,2, see (|21|). (j23j) and 
l|23|)-(|2^|l. 

From Lemmas HI [S] and El it follows that 

E(||m^)^0, (38) 

and 

E(\\Z«\\ Lix )<C. (39) 

Further on, from (|38j). respectively, (|39j) it follows that for any e > one 
can find a compact subset K e of and a bounded subset L e of x such 
that 



P [y n i K £ ) < e/2, P [Z n i L £ j < e/2. (40) 
On the other hand, Murat's lemma (see or Chapter 9 of |17j ) says that 

M £ := (K £ + L e ) n {X G : H^H^-i,- < C} 
is compact in H^. From (J3EJ), (EH) and (001) it follows that 

p(V^M e ) <£, 

uniformly in n, which proves the lemma. □ 
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5.6.2 Tartar's lemma and its consequence 

Lemma 8. Assume conditions fA\) and fflj) . Let Q n be a subsequence 
of the probability measures on y defined in \15\) . which converges weakly 
in the vague sense: Q™' =^ Q. Then Q is concentrated on the (vaguely 
closed) subset of Young measures satisfying That is, Q-a.s. for any 

two generalized Lax entropy/flux pairs (Si,F%) and (S 2 , F 2 ) and any test 
function if : [0, T] x T -> R, 

/ / <p(t,x)(v(t,x), SiF 2 - S 2 F 1 )dxdt= (41) 
Jo Jt 

f [ <p(t,x)((v(t,x), Si)(v(t,x), F 2 ) - (u(t,x), S 2 )(u(t,x), F x ))dxdt. 
Jo Jt 

Proof. First we prove (|41|) for twice continuously differentiable entropy/flux 
pairs. Due to separability of C([0, T] x T) it is sufficient to prove that for 
any two twice continuously differentiable Lax entropy/flux pairs (Si, F\) and 
(S 2 , F 2 ) and any test function ip : [0, T] x T -> R, (jUJ holds Q-almost-surely. 
So we fix (5*i, Fi), (6*2,^2) and ip. Note that 

:= dtSjir^xV+dxFjiT^x)) 
= d t {v n {t,x),S 3 )+d x {v n (t,x),F j ) 

J = 1,2. 

Due to Skorohod's representation theorem (see Theorem 1.8 of [3]) and 
LemmaElwe can realize the random Young measures is n (t,x;dv) and v(t,x;dv) 
jointly on an enlarged probablity space (H, A, P) so that P- almost- surely 

v n v, and {X™ : n' , j = 1,2} is relative compact in H^. 

So, applying Tartar's Div-Curl Lemma (see JH]) [20| . or Chapter 9 of [T7j ) 
we conclude that (in this realization) almost surely the factorization (|41j) 
holds true. 

Since P is compact, from Riemann's method of solving the linear hyper- 
bolic PDE ((7J) (see Chapter 4 of [Hj) it follows that generalized entropy/flux 
pairs are approximated pointwise by smooth ones. Thus the Tartar factor- 
ization ()41|) extends from smooth to generalized entropy/flux pairs. Hence 
the lemma. □ 
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The main consequence of Lemma |H1 is the following 

Proposition 3. Assume conditions FA\) and fflj) . Let Q n be a subsequence 
of the probability measures on y defined in \15)) . which converges weakly in 
the vague sense: Q n ' =^ Q. Then the probability measure Q is concentrated 
on a set of Dirac-type Young measures, that is Q(U) = 1. 

Proof. In view of Lemma |H1 this is a direct consequence of Lemma ^ □ 

Remark: This is the only point where we exploit the very special features 
of the PDE Q). Note that the proof of Lemma^relies on elementary explicit 
computations. In case of general 2x2 hyperbolic systems of conservation 
laws, instead of these explicit computations we should refer to DiPerna's 
arguments from [2] , possibly further complicated by the existence of singular 
(non- hyperbolic) points isolated at the boundary of the domain D. More 
general results will be presented in the forthcoming paper 

5.7 End of proof 

From Propositions |21 and |21 it follows that from any subsequence v! one can 
extract a sub-subsequence n" such that Q n " =^ Q and Q is concentrated 
on the set of Dirac-type measure valued entropy solutions of the Cauchy 
problem. From now on we denote simply by n this sub-subsequence. Refer- 
ring again to Skorohod's Representation Theorem we realize the Dirac-type 
random Young measures v^ x (dv) := 5-^^^^dv) and u tjX (dv) := S u ^ t ^(dv) 
jointly on an enlarged probability space (B,,A, P), so that v n — v almost 
surely and (t, x) \— ► u(t, x) is almost surely entropy solution of the Cauchy 
problem. From basic functional analytic considerations (see e.g. Chapter 9 
of pH) it follows that, in case that the limit Youg measure is also Dirac- 
type, the vague convergence v n — 1 v implies strong (i.e. norm) convergence 
of the underlying functions, 

f^u in L\ x . (42) 

So, we have realized jointly on the probability space (H, A, P) the empirical 
block average processes £ (t, x) and the random function u(t, x) so that the 
latter one is almost surely entropy solution of the Cauchy problem, and (|42|) 
almost surely holds true. This proves the theorem. □ 
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6 Appendix 

6.1 The logarithmic Sobolev inequality for random stirring 
of r colours on the linear graph {1,2,...,/} 

Let r > 2 be a fixed intger. For I € N we consider r-tuples of integers 
N = (iVi, . . . , N r ) such that 

N a > 0, a = 1, . . . , r and N x + ■ ■ ■ + N r = I, (43) 

I 

n l N := {w € {1, . . . , r} 1 : l{ Uj=a } = JV«, a = 1, . . . ,r}. 

Let 7rjy denote the uniform probability measure on Q l N : 

Aim) = — p , UL G ^4- 

The one dimensional marginals of 7r^y are 

1,1/ \ N a 

The random element of Q l N distributed according to tt 1 n will be denoted 
( = (Ci, (2, ■ ■ ■ , CO- Expectation with respect to 7rjy, respectively, Trjy will 
be denoted by Ejy (■■■), respectively, E^ 1 (•••)• Conditional expectation, 
given the first coordinate Ci will be denoted E^y( ■ • • Ki)- Note that 

E5v(/(C)|Ci = a) = E^i(/(a, C2, • • • , CO) 

where E^i( • • • ) stands for expectation with respect to (C2, • • • ,C0 dis- 
tributed according to n 1 ^} and, given N = (N\, . . . , N a , . . . , N r ) with N a > 
1, N a := (Ni,...,N a -l,...,N r ). 

Given a probability density h over (p, N ,ir l N ), its entropy is 

H l N (h) :=E l N (h(C) log h(()). 

Further on, for i,j G {1, . . . , 1} let Qij : fijy — > Q l N be the spin exchange 
operator 

Vj if fc = i, 
(®i,M) k = { if fc = j, 
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For / : Q l N — > K we define the Dirichlet form and the conditional Dirichlet 
form, given (j 

l-l 

D l N (f) ■■= -^E i 7V ((/(9 i , m C)-/(C)) 2 ), 
i=i 

l-i 

^(/|Ci) := -^E^((/(e M+1 c)-/(C)) 2 |Ci) 

■i=i 

= D i N c\(m,-)). 

The logarithmic Sobolev inequality is formulated in the following 

Proposition 4. There exist a finite constant N such that for any number of 
colours r, any block size I £ N, any distribution of colours N = (N\, . . . ,N r ) 
satisfying and any probability density h over (Q l N ,7r l N ), the following 

inequality holds: 

H l N (h) <Kl 2 D l N (Vh). (44) 

Remark: The proof follows (see also Due to exchangeability of 

the measures tt 1 n some steps are considearbly simpler than there. 

Proof. We shall prove the Proposition by induction on /. Denote 

H l N (h) 
W(l) :=supsup NK > . 

N h D L N {yfh) 

The following identity is straightforward 

H l N (h) = ^(^(higl^E^hxiO^ghxiOlCx)) 

+E# (e^MOICi) logEjv^ld)) , (45) 
where in the first term of the right hand side 

h m - h & 

hl ® '-W N (h(Q\CiY 

First we bound the first term on the right hand side of l|45jl . By the 
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induction hypothesis 

E^ 1 (E^(/t(C)|Ci) E^v (hi (C) log ^i(C)|Ci)) 

= E^ 1 (E l N (h(() | Ci)E l ~l (/>! (C) log hi (0) ) 



< 



w(l - i)e5v 1 (e5 v (/ 1 (C)|Ci)^ c 1 1 (V^i)) 



< W{l-l)D l N (Vh). (46) 

Next we turn to the second term on the right hand side of (|45|) . In order 
to simplify notation in the next argument we denote 

Qa ■= fe(i) := E l N (h(Ql {Cj=a} ). (47) 

It is straightforward that for any K < oo there exists a finite constant 
C = C{K) such that for any v E [0, K) 

vlogv <(v- 1) + C(\/v — l) 2 

and, furthermore, the constant C can be chosen so that for any v > K 

vlogv < Cv 3/2 . 

Hence, with the notation introduced in ()47|) . we get the following upper 
bound for the second term on the right hand side of (|45j) 

E^ 1 (E5 v (/ l (C)|C 1 )logE5 v (/ i (C)|Ci)) ^ea—^ — (48) 

Qa Qa 



a=l 
2 



a=l 



3/2 



We use the straightforward inequality 



^2 Qa ( ~ 1 ) 1 {2£«(1)<A'} - °- 



i L {2£Lii2<A'} 
fa / L ea — ' 

We choose K sufficiently large in order that Lemma 4.1 of [23 can be applied 
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to {1, 2, . . . , 1} 3 j i— > y/q a (J)/ 6a- Thus we obtain the upper bound 



3/2 



i 



Putting together and and returning to the explicit notation we 
obtain the following upper bound for the second term on the right hand side 
of (ESI): 



(e^v (^(C) I Ci) log E^v (^(C) |Ci) 



l-l r 



3=1 oc=l 



l-l r 



= c " l YH (V^C^iJ+iC)^^}) - \Jv l N (HC)hQ= a } 

j=l a =l 

= C"l £ Uv l N (h{Q h]+1 Q) - ^(MC))) 2 

< c"z Y, e n (y/HQn+iQ - = c"w l N (Vh) . (so) 

In the second step we used exchangeability of the canonical measures 71"^. 
In the last inequality we note that the map 

M + xM + 9(x i2/ )m(V^-^) 2 

is convex and we use Jensen's inequality. 

From (|l5|). (fl6|) and (|50|) eventually we obtain 

W(l) < W(l - 1) + C"Z, 
which yields (|Hj). □ 

6.2 An elementary probability lemma 

The contents of the present subsection, in paericular LemmaEland its Corol- 
lary n are borrowed form [23 • F° r their proofs see that paper. 
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Let (S7, 7r) be a finite probability space and iVi, i £ Z i.i.d. Sl-valued 
random variables with distribution ir. Further on let 

£ : n - R d , & : = ^(Wi), 

u : ->■ R, v i :=v(uji...,uji +m - 1 ). 

For cc G co(Ran(£)) denote 

r( , _ E 7r (t; 1 exp{E^ 1 A-^}}) 
L[X) - E7r (exp{A.^}}) m ' 

where co(-) stands for 'convex hull' and A G R d is chosen so that 

E^^expjA-^}}) 
E^(exp{A-| 1 }}) 

For IsNwe denote p/ain block averages by 

i=i 

Finally, let 6 : [0, 1] — > R be a fixed smooth function and denote 

M(b) := / 6(a) ds. 



x. 



Jo 

We also define the block averages weighted by b as 

I i 
< 6 ' t)i ■■= «>« : = yE^'/Ovi, 

j=0 j=0 
The following lemma relies on elementary probability arguments: 

Lemma 9. There exists a constant C < oo, depending only on m, on the 
joint distribution of (fj,£j) and on the function b, such that the following 
bounds hold uniformly in I G N and x G (Ran(£) + • • • + Ran(£))/Z: 
(%) IfM{b) = 0, then 

E(exp{ 7 V/(6, v)i} 1 1, =sc) <exp{C( 7 2 +7M)}- (51) 

fit; I/M(6) = 1 then 

E(exp{ 7 v^((6, t;),-T«6, 0/))} |?i = x) < exp{C( 7 2 + 7 M)}. (52) 

The proof of this lemma appears in \2'2\ . 
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Corollary 1. There exists a 70 > 0, depending only on m, on the joint 
distribution of (t>i,£j) and on the function b, such that the following bounds 
hold uniformly in I G N and x £ (Ran(£) + • • • + Ran(£))/Z: 
(i) IfMib) = 0, then 



x 



E^exp{ 7 oZ(&, v) 2 } 
(ii) IfMib) = 1 then 

E(e W { l0 l{(b,v) l -T((b,Oi)) 2 } 



< y/2. 



x 



< \/2. 



(53) 



(54) 



Proof. The bounds (|53|) and (|54"|) follow from (|5l|) . respectively, (|52|) by 
exponential Gaussian averaging. □ 

6.3 Proof of the a priori bounds (Proposition [^) 
6.3.1 Proof of the block replacement bound (j!8j) 

We note first that by simple numerical approximation (no probability bounds 
involved) 



/ v n {x) - T(T(x)) 2 dx - - \v n (j/n) - T(f (j/n)) 



< Cl 



n 2 a 



We apply Lemma |3] with 



V, 



v n (j/n)-T^ (j/n)) 



We use the bound 1)54(1 of Lemma ^ with the function b = a of (|13|) . Note 
that 7 = 70/ can be chosen in (jTTj) . This yields the bound (fl8|) . 

6.3.2 Proof of the gradient bound (fT9|) 

Again, we start with numerical approximation: 



1 n 

\d x v n (x)\ 2 dx--Y,\d x v n {j/n)[ 
We apply Lemma |3] with 



i=i 



V, = |^(j7n)| 2 . 
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We use now the bound (|53j) of Lemma Q with the function b = a' , where a 
is the weighting function from (|13|) . The same choice 7 = 70/ applies. This 
will yield the bound (fT9|) . 
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